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Abstract 

We reconsider QED radiative corrections (RC) to the it — > e + e~ decay width. One kind of RC 
investigated earlier has a renormalization group origin and can be associated with the final state 
interaction of electron and positron. It determines the distribution of lepton pair invariant masses 
in the whole kinematic region. The other type of RC has a double-logarithmic character and is 
related to almost on-mass-shell behavior of the lepton form factors. The total effect of RC for the 
7T° — ► e + e~ decay is estimated to be 3.2% and for the decay n — ► e + e~ is 4.3%. 

PACS numbers: 
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I. INTRODUCTION 



Rare decays of mesons serve as the low-energy test of the Standard Model. Accuracy of 
experiments has increased significantly in recent years. Theoretically, the main limitation 
comes from the large distance contributions of the strong sector of the Standard Model where 
the perturbative theory does not work. However, in some important cases the result can be 
essentially improved by relating these poorly known contributions to other experimentally 
known processes. The famous example is the Standard Model calculation of the anomalous 
magnetic moment of muon (g — 2)^ where the data of the processes e + e~ — > hadrons and 
t — > hadrons are essential to reduce the uncertainty. It turns out that this is also the case 
for the rare neutral pion decay into an electron-positron pair measured recently by the KTeV 
collaboration [l] and reconsidered theoretically in 

The measured branching is [l| 

^KTeV _^ e + e - jXD > 0.95) = (6.44 ± 0.25 ± 0.22) ■ 10" 8 , (1) 

where the kinematic cut over the Dalitz variable xd = {p+ + P-) 2 /M 2 , v 2 = Am 2 /M 2 < 
xd < 1? was used in order to suppress the Dalitz decay events 7r° — > e + e~ r y. Then, the 
important step in extraction of the branching consists in correct treatin g th e radiative cor- 



rections (RC) to the process which has been considered earlier in [3[ and [llj. Extrapolating 



the full radiative tail beyond xn > 0.95 and scaling the result back up by the overall RC 
leads to the final result lj 

^ TeV (vr -> e + e-) = (7.49 ± 0.29 ± 0.25) • 10" 8 , (2) 



where the leading order radiative corrections have been taken into account [3j|. It is the 
motivation of our paper to revise the calculation of QED RC to the 7r — > e + e~ decay width. 

In the lowest order of QED perturbation theory (PT), the photonless decay of the neutral 
pion, 

2 j\ 4-2 2 _2 



7T (g) -> e (p_) + e + (p + ), q 2 = M , p\ 



m 



(M meson mass, m lepton mass) is described by the one-loop Feynman amplitude (Fig. la) 
corresponding to the conversion of the pion through two virtual photons into an electron- 



positron pair. The normalized branching ratio is given by 



Bn (n° — > e + e ) n / am \ 2 . . , , .„ 
- C e ' ] = Bin ^ -,-,) = 213 (rf?) ' M I • < 3) 



B (tt° — > 77) 



FIG. 1: Set of the lowest order QED RC to tt° — > e + e process: virtual corrections. 



where (3 = a/1 — z/ 2 , 5 (7r° — > 77) = 0.988 and the reduced amplitude is 
d 4 k {qkf - q 2 k 2 



2 

M 2 



F,{-k\-{q-kY), (4) 



m 2 (k 2 + it) [(q — k) 2 + ie] [(p_ — A;) 2 — m 2 + z'e] 
with the pion transition form factor F n (—k 2 ,—q 2 ) being normalized as ^(0,0) = 1. The 
imaginary part of A (q 2 ) can be found in a model independent way [5] 

'l + P 
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\mA{M z ) = In 



(5) 



2/3 \l-P, 

while the real part is reconstructed by using the dispersion approach up to a subtraction 
constant 



ReA(M 2 ) = A(0) + - 



h - In 2 

12 4 



P 



1-/5 

Usually this constant, containing the nontrivial dynamics o: 



(6) 



the process, is calculated 



tne pro 



within different models describing the form factor ^(/c 2 ,^ 2 ) [2j, [6|, [7|, [8j. However, it has 
recently been shown in [2| that this constant may be expressed in terms of the inverse moment 
of the pion transition form factor given in symmetric kinematics of spacelike photons 



^( g 2 = 0)=31n(^ 



3 
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(t,t) 



t 



F * . 
dt- 1 7 



(M) 



t 
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Here, \i is an arbitrary (factorization) scale. One has to note that the logarithmic dependence 
of the first term on \x is compensated by the scale dependence of the integrals in the brackets. 
The accuracy of these calculations are determined by omitted small power corrections of the 
order O(^p-) and O(^aL) in the r.h.s. (jSJ), where A < M p is the characteristic scale of the 
form factor -F 7r7 * 7 » (t, t) and L is the large logarithm parameter 

'M 2 \ . 



In 



In 



P 



For the decay n° — > e + e _ one has L pa 11.2. 

By using the representation (JTj), and the CELLO [9j and CLEO [lOj] data on the pion 
transition form factor F^^? (t, 0) given in asymmetric kinematics the lower bound on 
the decay branching ratio was found in This lower bound follows from the property: 
F nl * y * (t, t) < F 7r7 . 7 . (t, 0) for t > 0. It considerably improves the so-called unitary bound 
obtained from the property \A\ 2 > (Im.,4) 2 . Further restrictions follow from QCD and allow 
one to make a model independent prediction for the branching 

^Thcor ^0 ^ e+e -) = ±QA y ^ (g) 

which is 3.3a below the KTeV result (J2J). The main source of the error in (jHJ) is defined 
by indefiniteness in the knowledge of the pion form factor i*Vy* 7 * (t,t) The discrepancy 
between OH]) and (j2j) requires further attention from experiment and theory to this process 
because there are not many places where experiment is in conflict with the Standard Model. 

Considering the higher orders of QED PT, there are two sources of RC to the width of 
the 7T° — > e + e~ decay (Fig. 1). One of them has a renormalization group origin and can be 
associated with the final state interaction of electron and positron. The relevant contribution 
corresponds to taking into account the charged particle interaction at large distances. The 
other is of double-logarithmic character and is related to short distance contributions. Let us 
note that the branching ratio ([3]) is proportional to the electron mass squared and thus the 
Kinoshita-Lee-Nauenberg theorem of cancellation of mass singularities in the limit m — > 
is not violated. 



e (p_) 






FIG. 2: RC due to soft photon emission. 



The first kind of corrections was considered in |3j]. Later, the effect of the higher order 



RC was estimated in [HI] by using the exponentiation of soft photon contributions, which 
is essentially equivalent to the Yennie, Frautchi-Suura factorization procedure. Considering 
only the two-virtual-photon conversion to a lepton pair, it is originated from the box type 
Feynman amplitude (Fig. le) as well as the contribution from the emission of real photons 



by leptons (Fig. 2) and produces the single-logarithmic enhanced terms (~ L) which are 
described by the lepton nonsinglet structure function method [ijj]. It was shown in jll] that 
the soft photon emission can drastically change the results obtained at the Born level when 
the invariant mass of leptons is close to the pion mass. 

Furthermore, we find an additional source of RC which is of the so-called "double- 
logarithmic" (DL) nature (aL 2 /n ~ 1). This kind of asymptotics was intensively inves- 



tigated in the 70s in a series of QED processes 
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14j. The DL type contribution to the 



decay width was not considered earlier for the 7r — > e + e decay. 



II. LARGE LOGARITHM REGIME AND DOUBLE-LOGARITHMIC CORREC- 
TION 

First, it is instructive to reproduce the results discussed above in a simple and physical 
way. For this aim we note that the main contribution to the real part of A (M 2 ) comes from 
the kinematic region of loop momenta corresponding to the intermediate virtual electron (or 
positron) close to the mass shell (Fig. la). Really, by changing the integration variable in 
(1) as qi = k = p_ — k, q 2 = q — k = p + + K and omitting terms of order 0(k 2 /M 2 ) we 
can rewrite the amplitude in the Born approximation Ao (M 2 ) as 

1 fd A K M 2 



A (M 2 ) 

2 J in 2 (k 2 — m 2 + ie)((K — p^) 2 + ie)((n + p + ) 2 + ie) 

Let us find the real part of the amplitude within the leading logarithmic accuracy that 
corresponds to the restriction of the kinematic region by conditions 



2 I 1 2 

m ~ \k\ 



«(| gi 2 |,|g 2 2 |)«M 2 . (9) 



To this end one performs the substitutions 
d 4 n — in~K 2 d\lt \dO 



A + ie 2co 

q\ = -2Mk u, q\ = 2Mk (u + U cos 9) , 



[eW + QNlb, ( 10 ) 



/=^- o A 2 2 a L m 2 1-&;COS0 

u=VK 2 + m 2 , A = k -m, fj w = \ I u = 

V uo z 2 

where 9 is the angle between the directions of electron momentum (the rest frame of the 
initial pion implied) and the 3-momentum of the virtual electron. Let us note that by 
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kinematical reasons in the region of maximal contribution the signs of q\ and q\ must be 
opposite. Performing the angular integration we obtain the leading term of ([6]) 



ReA M 2 = / — ^ = / ^— In — « -L . 11 

y 7T J ^g 2 g 2 Jm ^ m 2 4 

Then, let us consider the vertex type RC (Fig. 1 c,d). The lowest order evaluation arising 
from the diagrams of Figs, lb) and lc) leads to the correction T(qi, n,p_) = 1 — |^iy(g 2 , n 2 ) 
with 

r i 2 2n , kill 1,2 3 , \Ql\ , ^ , 1 , 171 a ( 2\3vT 2 

Mgi,-) = ln^ln M + -ln R - - In — + y + - - In y - 9 — , 

which is consistent with the result of similar calculations in The last term arises from 
a renormalization procedure. A similar contribution Iv(q%, ft 2 ) comes from the diagrams of 
Figs, lb) and Id). 

Let us consider now the box-type diagram (see Fig. le) and demonstrate the calculations 
in more detail. The corresponding contribution to the amplitude has the form 

^Re [*h a. 

2 J i-K 2 (kf — A 2 + it) {{p~ + k{) 2 — m 2 + ie) ((p + — ki) 2 — m 2 + it) ((k + ki) 2 — m 2 + ie) 

(12) 

where 

N = u(p_)7 M (p_ + ki + m)7 A (« + h + m)^(-p + + k x + m)j /J ,v(p + ). 
In the leading kinematic region, where k\ and k are small, we can reduce the numerator to 

N « —2mM 2 u{p_)^xl'yv{p + ). 

The calculation of the scalar 4-denominator integral is standard: by using the Feynman 
parametrization and performing loop momentum integration we arrive at 

f 1 f l f 1 z 2 dz 
—2mM 2 A dx ydy 

Jo Jo Jo 



(Az 2 + Bz + C) 2 ' 



with 



A = (y Px - yn) 2 ; B = -yA - A 2 ; C = A 2 , (13) 
p x = xp + — xp_, pi = m 2 — M 2 xx, x = l — x, y — 1 — y. 
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First, we perform the integration in z 

z 2 dz 2C + B 



2C + B + ^/R 

(Az 2 + Bz + C) 2 ~ (A + B + C)R~Rl n 2C + B- Vr' 

R = B 2 - AAC > 0. 



(14) 



Calculating the y-integral of (fT4j) results in 



1 A 2 1 p 2 
-In— + ln^-ln 

2 \ z m z 2 m z 



q\x + q\x 



Integration in x by using 
" 1 dx 



Re 



l p 2 x + ie M 2 
leads to the correction t^/r with 

Z7T 



2 L, R e r-^-ln^ + " 







A,(L 2 --7T 2 ) 

M 2V 3 ; 



Jfl (g?, g 2 2 ) = --L 2 - 2(L - 1) In - - L(L 1 + L 2 ) -Ul x - L 2 ) 2 + -tt 2 



(15) 



where 



L12 = In 



1 01,2 1 



m 



2 ' 



Finally, one needs to integrate over photon momenta q\ and q 2 - Again, the logarithmically 
enhanced contribution comes from the kinematic regions m < uj < M/2, cos 6 — > ±1 (see 
definitions in (fTUl) ). The Born amplitude (one- loop) and the lowest-order radiative correction 
to it can be written as (we take into account the equal contributions of regions \qf\ <C |g|| 
and \q 2 \ < \qf\) 

du 



UJ J m2 U 



1 + ^(Iv (ql m 2 ) + l v (q 2 2 , m 2 ) + I B (qj, q 2 2 )) 



(16) 



with (we put here A « m 2 



I v (ql m 2 ) + I V (ql m 2 ) + I B (ql q 2 2 ) = -\ L 2 - hi - \l 2 - 2(L - 1) In ^ + \l + \l 



- l 2 u - l-L + - (I - 1) \l u + — + 1 



2 2 

where we use substitutions ffTOl) and introduce the notation 



7T 
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(17) 



I = In 



L = In 



1 — Puj cos 9 



(18) 
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Integration of (Tl6|) leads to (we keep only terms of order L 2 , L and LP) 

Rvirt 



R{) 

&virt 



1 ~\~ O v i r t) 
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, , , m 3 7i 

— L — 2(L — 1) In h -L H 

24 V A 4 6 



(19) 



Consider now the real photon emission corrections. One can distinguish two mechanisms 
of the radiative decay tt° — > e + e~ r y. One of them, the so-called Dalitz process, corresponds 
to decay mode of the pion to real and virtual photons with a subsequent decay of the virtual 
photon to the e + e~ pair. The corresponding contribution to the width is not suppresse d by 
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17| 



lepton mass and provides an important background to the tt — > e + e process |3 
However, the Dalitz matrix element squared and its interference with the double virtual 

3 2 n 

photon amplitude are suppressed by ~ (1 — xd) and ~ (1 — xd) as — > 1 [3f (Fig. 3). 
This results in a negligible (of order 0.02%) interference contribution integrated in the region 



19j. 



of interest for this measurement, 0.9 < xd < 1 

Another mechanism consists in creation of a lepton pair by two virtual photons with 
emission of real photon by a pair components. For emission of a soft photon (with energy 



to not exceeding Ae << 4f in the pion rest frame) the standard calculations [l3j give 



Jsoft 



a 

7Y 



2Ae 



2(L-l)ln — + 2(L-l)lny + -L 
Emission of a hard photon was investigated in with the result 



Ti" 

~3~ 



(20) 
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hard 



a 
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(21) 



In order to find the distribution over a lepton pair invariant mass, the adequate way is to 
use the method of structure functions [l2|] based on the application of the renormalization 
group approach to QED. Here, the nonsinglet structure function is associated with a final 
state fermion line. In partonic language it describes the probability for a fermion to stay a 
fermion. Omitting the events with creation of more than one lepton, the nonsinglet structure 
function of electron relevant to the process valid at all orders in perturbation theory is 1 



F(x D ) =b(l- xjo)"- 1 (l + h^j - + x D ) + O (b 2 ) 



(22) 



1 Note that the value of the i^-factor in (|22| , K = 1 + 1 — , differs from that obtained in [12| for the process 
of hadron production in single-photon e + e~ annihilation channel. 



where b = — (L — 1). To the lowest order in b and in the region xd » v 2 = 4m 2 /M 2 , the 
above expression is in agreement with the leading order expression j^] 

1 dR^S m (x D ) a 1 j\ , 2 , , (l + f3 x 



Ro d,XD 7T 1 — Id 



|(l+4)ln(^)-2x^}, (23) 



where 



Thus we arrive at the differential rate 

1 dR* c (x D ) 



u 2 



X D 



JF(x D ), (24) 



Ro dxo 

where the normalization factor J takes into account to total RC. The distribution (1241) . 
shown in Fig. 3, in contrast to (1231) is free of nonintegrable singularity at xd = 1- So, we 
see the importance of taking into account the higher orders of perturbation theory. 
Adding the virtual and real photon emission contributions we finally obtain 



Rn 



J / F(x D , L)dx D = 1 + a virt + a so ft + awd 
Jo 



1-* 

7T 



' 1 r2 3 r vr 2 21 
24 L + 4 L "T + T 



0.968. (25) 



The total radiative correction contains the large logarithm term ~ L because the 75 
current is not conserved. Compared with [3J we provide a more detailed analysis of the 
effective tt° — > e + e~ vertex revealing its DL structure 2 . However, numerically the ratio of 
the total RC corrections to the lowest level rate estimated in §25§ is very close to —3.4% 
found in [3] and used by the KTeV Collaboration in their analysis. 

The branching ratio of the pion decay into an electron-positron pair has been measured 
by the KTeV collaboration in the restricted kinematic region in order to avoid a large 
background from the Dalitz process dominating at lower values x&. By using the distribution 
( 1221 we can estimate the factor of extrapolation of the full radiative tail beyond xd > 0.95 
as /0.95 = 1.114. With this factor and scaling the result back up by the overall radiative 
correction ([23]) we confirm the result d2J) obtained by the KTeV Collaboration. 



2 It is naturally to expect the existence of DL-type RC in higher orders of perturbative theory. We do not 
touch this problem here. 
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FIG. 3: Distribution of different contributions over x = -jjt- solid line - Dalitz mechanism con- 



tribution, dashed line - inner bremsstrahlung contribution (|24p . dotted line - interference of these 
two mechanisms. 

III. CONCLUSIONS 

In this work, we reconsidered the contribution of QED radiative corrections to the tt° — > 
e + e~ decay which must be taken into account when comparing the theoretical prediction (jSJ) 
with experimental result (JTJ). Comparing with earlier calculations ^j], the main progress is 
in detailed consideration of the 7*7* — > e + e _ subprocess and revealing of dynamics of large 
and small distances. The large distance subprocess associated with final state interaction 
produces the terms linear in the large logarithm parameter L. The double logarithmic 
contributions (~ L 2 ) correspond to configurations when the particles in the loop are highly 
virtual. The total result is in reduction of the normalization factor by 1 — J ~ 0.032. 
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Occasionally, this number agrees well with the earlier prediction based on calculations [3|] 
and thus we confirm the KTeV analysis of RC factors. So our main conclusion is that 
taking into account of radiative corrections is unable to reduce the discrepancy between the 
theoretical prediction for the decay rate ([H]) and experimental result (TS]). Further independent 
experimental and theoretical efforts are necessary. Note, that if the discrepancy will stand 
as it is, then the possible explanation of the effect is due the contribution to the decay width 



of low mass (~ 10 MeV) vector boson appearing in some models of dark matter 211 ] . 

The question about the corrections is also important for the decays of r\ and K mesons 
to fi + fi~ and must be taken into account for the analysis of experimental data. However, for 
these decays the large logarithm parameter does not arise. The analysis of the lowest order 



RC to the decay width of kaon to rmion pair was given in 
of 



181 ] . Unfortunately, the result 



181 ] explicitly depends on infrared singularities and cannot be used in practice. Thus, 
in order to extract the information about P — > 11 decays (where P = Kl, 7To, 77, . . . and 
I = e,fi) calculated within the frame of definite models, the RC must be taken into account 



when working with experimental data, e.g., for the process 77 — > 20J]. For the process 

77 — > e + e~ our results are applicable and we estimate RC at the level 4.3%. 

Finally, we have to emphasize that the role of RC is rather timely because of the growing 
accuracy of modern experiments. 
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